Dynamics of lane formation in driven binary complex plasmas 
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The dynamical onset of lane formation is studied in experiments with binary complex plasmas 
under microgravity conditions. Small microparticles are driven and penetrate into a cloud of big 
particles, revealing a strong tendency towards lane formation. The observed time-resolved lane 
formation process is in good agreement with computer simulations of a binary Yukawa model with 
Langevin dynamics. The laning is quantified in terms of the anisotropic scaling index, leading to a 
universal order parameter for driven systems. 
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The formation of lanes is a ubiquitous phenomenon 
occurring in nature when two species of particles are 
driven against each other. When the driving forces are 
strong enough, like-driven particles form "stream lines" 
and move collectively in lanes. Typically, the lanes ex- 
hibit a considerably anisotropic structural order accom- 
panied by an enhancement of their mobility. The phe- 
nomenon is most commonly known from pedestrian dy- 
namics in highly populated pedestrian zones [l| , but also 
occurs in very different model systems of driven particles, 
such as colloidal dispersions 0] , lattice gases [5j] and 
molecular ions Q. Lane formation is a genuine nonequi- 
librium transition ^ which depends on the details of the 
particle interactions and their dynamics 

Recently, the particle laning was also observed in com- 
plex plasmas [sj. In fact, complex plasmas [^[l^ provide 
a very important intermediate dynamical regime that is 
between classic undamped fluids and fully damped col- 
loidal suspensions: In complex plasmas, the "atomistic" 
dynamics associated with the interparticle interaction 
is virtually undamped whereas the large-scale hydrody- 
namics is determined by friction. 

In this Letter we report on comprehensive experimen- 
tal studies of lane formation in complex plasmas, that 
were carried out under microgravity conditions with the 
radio frequency (rf) discharge chamber PK-3 Plus 
The motivation for this research is threefold: First, we 
demonstrate that complex plasmas are indeed an ideal 
model system to study nonequilibrium phase transitions 
such as laning. Second, the experiments enable us to 
investigate the dynamical onset of lane formation in de- 
tail. Third, we achieve a quantitative understanding of 
the structural correlation during the onset of laning by 
comparison with particle-resolved Langevin simulations. 
Based on the anisotropic scaling index analysis of the ob- 
tained data, we suggest a universal order parameter for 
non-equilibrium phase transitions in driven systems. 

Experiments. A series of dedicated experiments was 
carried out in the PK-3 Plus rf discharge chamber on the 
International Space Station. Microgravity is a crucial 
factor which allows us to form stable spheroidal complex 



plasma clouds. When small particles are injected into 
a cloud of big ones, the force field pulls small particles 
towards the center, thus making such systems perfectly 
suited to study lane formation. Our experiments were 
performed for various combinations of "big" and "small" 
monodisperse particles (2.55, 6.8, 9.2, and 14.9 /im diam- 
eter for "big", and 1.55, 2.55, and 3.4 ^m for "small"), 
with different neutral gases and pressures (argon between 
10 and 60 Pa and neon at 60 Pa, to control the friction 
rate) , and for different rf discharge powers (to control the 
screening length). 

Figure [1] shows a characteristic example of lane forma- 
tion observed in experiments with 3.4 /im and 9.2 ^va par- 
ticles at pressure of 30 Pa. When a fraction of individual 
small particles enters the interface of the fairly homoge- 
neous cloud formed by big particles, the subsequent pene- 
tration is accompanied by a remarkable self-organization 
sequence: (a) Big particles are pushed collectively by the 
inflowing cloud of small particles, the latter form strings 
drifting on average along the force field, (b) As the par- 
ticles approach the center of chamber, the field decreases 
and the strings organize themselves into larger streams. 
At the later stage, when the field almost vanishes, the 
streams merge to form a spheroidal droplet with well- 
defined surface, indicating the transition to the regime 
when surface tension plays the primary role. In this Let- 
ter we are specifically interested in the first two stages 
[T^ . It is noteworthy that during the stage (b) big parti- 
cles also form well defined strings. Small and big particles 
create an "array" of interpenetrating strings, where small 
particles are drifting to the right and big ones are almost 
immobile. After the flux of small particles is exhausted, 
the big-particle strings slowly dissolve [13]. 

Computer simulations. First, the distribution of the 
characteristic parameters of the discharge plasma, such 
as plasma density, electron temperature, and electric 
fields, was deduced from 2D simulations of the PK-3 Plus 
discharge chamber. We used the plasma fluid code (de- 
scribed in Ref. Jjj), which provides a self-consistent cou- 
pling of dust species to the discharge plasma, including 
the particle charging, plasma absorption, redistribution 
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FIG. 1: Lane formation in complex plasmas. A short burst 
of small (3.4 /im) particles is injected into a cloud of big 
(9.2 nm) background particles (close to the midplane of the 
chamber, indicated by horizontal dashed line) . Small particles 
are driven towards the center, stages of (a) initial lane forma- 
tion and (b) merging of lanes into larger streams are shown. 
Particles are illuminated by a thin laser sheet of ~ 0.35 mm; 
each figure is a superposition of two consecutive color-coded 
images (l/50th s apart, red to blue), the time difference be- 
tween them is ~ 1.2 s. At the stage (b) big particles also 
form well defined lanes that can be identified as strings of red 
dots. The frame indicates the region used for the analysis of 
big-particle dynamics. 



of volume charges/ambipolar fields, etc. These simula- 
tions suggest that in the midplane of the chamber the 
cloud of big particles is "self-confined" at its edges due 
to the self-consistent plasma field, as illustrated in Fig. 13 
of Ref . [3] , whereas inside the cloud the electric and ion 
drag forces practically compensate each other. On the 
other hand, since these two forces have different scaling 
on the particle size, the net force on small particles is 
nonzero inside the cloud. For the experiment shown in 
Fig. [U the simulations yield a force of fs = 0.3 ± 0.1 pN 
pushing a small particle towards the center. 

Next, particle-resolved molecular dynamics (MD) sim- 
ulations were performed on the Langevin level [15| for 
a binary mixture of 5759 small and 12287 big particles. 
The simulation box with periodic boundary conditions 
has dimensions of 4.4 cm in x-direction and 0.8 cm in y- 
and z-dircctions. The particles interact via a Yukawa pair 
potential with a screening length A = 100 ^m (based on 
results of our plasma discharge simulations) and effective 
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FIG. 2: Lane formation in MD simulation corresponding to 
the experiment shown in Fig. [T] Two snapshots illustrate (a) 
the initial injection stage and (b) the steady state. 



charges Zs = 3000 e (based on experiment and Zb = 
8117 e, proportional to the respective particle diameter, 
CTs = 3.4 ^m and ah = 9.2 fim. The mass density of 
the particles is 1.5 g/cm'^ and the corresponding friction 
rates are Vs = 250 s~^ and = (o's/crb)^'s = 92.4 s~^. 
The mean interparticle distances are deduced from the 
experiment, Ag — 464 iim (before the penetration) and 
Ab = 493 fim, the temperature is T = 0.024 eV. 

Our plasma simulations show that big particles do not 
experience an external force in the bulk. Therefore, in 
our MD simulation we confine them in x-direction by a 
parabolic external potential at two edges (with a width of 
2.2 cm) and adjust the confinement strength, so that the 
measured interparticle distance Ab is reproduced. Simi- 
larly, a portion of small particles, separated from the big 
particles, was prepared (with a width of 0.7 cm). Then 
the constant driving force fs was instantaneously applied, 
leading to penetration of small particles into the cloud of 
big ones. Simulation snapshots are shown in Fig. ^ re- 
vealing a qualitative agreement with the experiment. 

Anisotropic scaling index and order parameter. In or- 
der to identify and quantify the string-like structures in 
our experimental data and the simulations, a suitable or- 
der parameter has to be employed that is sensitive to the 
changing particle structures. Conventional approaches, 
e.g., binary correlation or bond orientation functions, 
Legendre polynomials, etc. turned out to be too insensi- 
tive. Much more satisfactory results were obtained by im- 
plementing an anisotropic scaling index method - a local 
nonlinear measure for structure characterization. This 
method has already been used to characterize electrorhe- 
ological cornplex plasmas 17|, large-scale distribution of 
or bone structure [l9| . 

1,...7V, 



galaxies [15 

For a given set of particle positions, {r;}, i 
we define a local density p{ri,R) — J2f=i s{dij/R), where 



dij — 



Ti — Tj I and s is a certain shaping function charac- 



terized by the spatial scale R. The scaling index a is the 
logarithmic derivative of the density with respect to the 
spatial scale, a = d\ogp{ri,R)/d\ogR. Hence, a{ri,R) 
characterizes the dimensionality of the local structure 
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around point r^, in the vicinity determined by the scale 
R. For example, a{ri,R) ~ 1 means that the local struc- 
ture is close to a straight line at the spatial scale R, for 
a{ri, i?) ~ 2 it is an element of a plane, and so on. Using 
the Gaussian shaping function s = e"'^'^'^/^^ we derive 
for the scaling index, 



a{ri,R) = 



1^3 = 1 



(1) 



Thus, the spatial scale R is an important mesoscopic 
measure of the local environment. In the range of rel- 
evant scales a is practically independent of R. 

In order to characterize anisotropic structures, we use 
a "stretch metric" for the distance measure dij. On a 
2D plane, the metric is determined by the aspect ratio 
e (> 1), which is the relative stretching of two principal 
axes, and by the unit vector u = (cos 0, sin 0) in the di- 
rection of stretching. Then the resulting anisotropic scal- 
ing index, a{ri, R,9), can be directly obtained from Eq. 
([!]) . We propose a "uniaxial vector characterization" : 
Each point is associated with the unit vector which 
points to a "preferred" direction of the local anisotropy. 
This direction is determined by the angle 9i at which the 
"anisotropic contrast" a(ri, i?, 9i + tt/2) — a{ri,R, 9i) is 
maximized. The directions and — are equivalent, so 
that below they are defined for the range < 9i < ^. 

Thus, each point can now be considered as a uniaxial 
"molecule" (simple rod) with the direction u, . Therefore, 
the global laning on a 2D plane can be characterized with 
the second-rank tensor Qap = 2N~^ 'l2iLi Ui (X> — Safi, 
analogous to that used to quantify order of the nematic 
phase. The direction of the global laning, (u), is then 
the eigenvector ( "nematic director" ) corresponding to the 
largest eigenvalue of QafS, which in turn is the laning or- 
der parameter, S. Obviously, S = 1 for a perfect align- 
ment and S = for a disordered phase, when individual 
vectors Ui are uncorrelated. We finally define the global 
laning angle O via the relation cos 8 — (u) • e^;. 

Comparison of experiments and simulations. For the 
analysis of the MD simulations, we divided the volume of 
the simulation box in z-direction (perpendicular to the 
driving force in x-direction) into several slabs of 0.35 mm 
width, which is about the thickness of the laser sheet 
used to record the experimental data. Hence, the aver- 
age number of particles and their average density, as well 
as the magnitude of their fluctuations in this "reduced" 
simulation data set were similar to those in the experi- 
ment. The obtained results were analyzed and compared 
using the anisotropic scaling index method. The relevant 
range of spatial scales R was about (2.5 — 4)Ab for big 
particles and (3 — 5)Ab for small particles, the anisotropic 
aspect ratio e was 5 and 7, respectively. 

Discrimination of big and small particles in the ex- 
perimental data was performed in terms of their veloci- 
ties: During the stages of the lane evolution illustrated 
in Fig. [U small particles drift relatively fast with respect 
to big particles towards the center of the chamber. They 



gradually slow down, until at the later stage their veloci- 
ties become too low and then small and big particles are 
no longer distinguishable. 

(i) Small particles. We first identified boundaries of 
the self-contained big- and small-particle clouds. Over- 
lap of the two clouds was a natural choice for the region 
to analyze the dynamics and structural evolution of small 
particles, as they penetrate the cloud of big ones. There 
are two distinct phases characterizing formation and evo- 
lution of small-particle lanes [approximately correspond 
to Figs[T](a) and[ljb), respectively]: An "injection stage" 
(I) , starts at the moment when first small particles pene- 
trate the cloud of big ones. During this phase, the num- 
ber of small particles for the analysis increases from zero 
to the average number per slab, and therefore significant 
fluctuations are possible due to poor statistics. After 
about one second, there is a crossover to a "steady-state 
stage" (II), when the average number of small particles 
remains constant and the driving force can be consid- 
ered constant as well. In the simulations the duration of 
steady-state phase is sufficiently long (two seconds), due 
to properly chosen length of the simulation box. In the 
experiment, however, this phase is 2-3 times shorter, be- 
cause the driving force and hence the particle velocities 
decrease as they approach the center of chamber. 

The order parameter Ss{t) calculated for small parti- 
cles is plotted in Fig. . Although in the beginning of 
stage I it exhibits significant fluctuations due to poorer 
statistics, one can see that the formation of small-particle 
lanes is practically "instantaneous" at timescales corre- 
sponding to the experimental frame rate (50 frames/s). 
The magnitude of the order parameter in the experiment 
is almost twice as low as that in the MD simulation. We 
believe that this discrepancy is due to the fact that the 
discrimination procedure allows us to identify only 70- 
80 % of small particles in the experimental data, which 
results in the artificial "thinning" of the small-particle 
lanes. The random elimination of 30 % small particles in 
the MD data decreases Ss down to the experimental level, 
clearly supporting this hypothesis. Figure[3Jb) shows the 
evolution of the global laning angle Qsit) for small parti- 
cles, which exhibits narrow dispersion and demonstrates 
that the "nematic director" practically coincides with the 
driving vector. Note the increasing broadening and de- 
flection of 8s (t) seen in experiment at the "steady-state" 
stage II [coinciding with slight decrease in <S's(t)], which is 
related to convergence of small-particle lanes towards the 
chamber center as the driving force gradually decreases. 

(ii) Big particles. In order to diminish possible influ- 
ence of the boundary effects, we defined fixed regions in 
the bulk of the big-particle cloud, as indicated in Figs [1] 
and [21 In terms of the size and shape, these regions ap- 
proximately correspond to the "overlap regions" used for 
the analysis of small-particle dynamics. 

The evolution of the order parameter and of the global 
laning angle for big particles is shown in Figs and 
[31[d), respectively. Initially, there is no anisotropy in the 
simulation, whereas in the experiment — 0.05 and 
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FIG. 3: Dynamics of laning. Shown are the evolution of 
the "nematic" order parameter for small and big particles, 
Ss (a) and Sh (c), respectively, as well as the corresponding 
global laning angle, Os (b) and 0b (d), as obtained from the 
anisotropic scaling index analysis of the experiment (crosses) 
and MD simulation (dots). For small particles, injection stage 
I and steady-state stage II are indicated. 



Ob ^ 0, due to a weak inhomogeneity in the big-particle 
density. Once a cloud of small particles reaches the fixed 
region used for the analysis, starts growing and the an- 
gular distribution narrows around Ob — 0, due to the in- 
creasing number of small particles causing the formation 
of big-particle strings. Then Sb reaches a maximum and 
starts falling off, reflecting the onset of string relaxation 
when small particles leave the region. The initial relax- 
ation occurs at a characteristic timescale of ^ 1 s which 
is an order of magnitude shorter than the self-diffusion 
timescale for the big particles (~ m^VbA'^/T). Note that 
after this rapid relaxation Sh (t) tends to some intermedi- 
ate plateau, and laning angle Qh{t) keeps the anisotropy, 
indicating that the structural relaxation is apparently not 
complete. This suggests that the ultimate equilibration 
might involve some metastable states. 

Thus, binary complex plasmas provide us with cru- 
cial new insights into very important dynamical regime 
of laning, mediating classic undamped fluids and fully 
damped colloidal suspensions. By combining the experi- 
mental studies and the particle-resolved Langevin simu- 
lations, we investigated the dynamical onset of lane for- 
mation in driven complex plasmas. Furthermore, we pro- 
posed a universal order parameter for the characteriza- 
tion of non-equilibrium phase transitions in driven sys- 
tems. The approach is based on the anisotropic scaling 
index analysis that is exceptionally sensitive to symmetry 
changes occurring in particle ensembles. The use of such 
order parameter could be invaluable in studying the onset 
of non-equilibrium phase transitions. In particular, this 
might help us to shed the light into the principal issue 
about the order of phase transition, characterize possible 
universality, identify dynamical regimes of structural re- 
laxation, etc. As the immediate next steps, one can think 
of employing the proposed approach to investigate laning 
in periodically driven [l^l and crystalline 21] systems. 
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